The role of E1-E2 interplay in multiphonon Coulomb excitation 
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In this work we study the problem of a charged particle, bound in a harmonic-oscillator potential, 
being excited by the Coulomb field from a fast charged projectile. Based on a classical solution to 
the problem and using the squeezed-state formalism we are able to treat exactly both dipole and 
quadrupole Coulomb field components. Addressing various transition amplitudes and processes of 
multiphonon excitation we study different aspects resulting from the interplay between El and E2 
fields, ranging from classical dynamic polarization effects to questions of quantum interference. We 
compare exact calculations with approximate methods. Results of this work and the formalism we 
present can be useful in studies of nuclear reaction physics and in atomic stopping theory. 
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I. INTRODUCTION 

Recent studies of Coulomb-induced breakup of weakly 
bound nuclei indicate a considerable reduction of the dis- 
sociation probability in comparison with the prediction of 
the first-order Born calculation p| . The reduction is asso- 
ciated with dynamic polarization, where the quadrupole 
field creates a polarization that influences the effects in- 
duced by the dominating electric dipole. This is a well- 
known phenomenon in atomic stopping theory, where the 
stopping power for charged particles deviates from the Z 2 
dependence of Bethe's formula @], Z being the charge of 
a particle. Measurements, as for example, of the stop- 
ping power for protons and anti-protons |3|], indicate the 
presence of a Z 3 correction. Over the years a substantial 
progress has been made in measuring and understanding 
of this phenomenon 0, H IS II II S Il3 • 

In this work we consider the Coulomb excitation of 
a charged particle bound by a harmonic-oscillator po- 
tential. The time-dependent electric dipole (El) and 
quadrupole (E2) fields of a projectile excite the parti- 
cle, and within this approximation the quantum prob- 
lem is treated exactly. This problem has been studied in 
the past, starting from the classical treatment by Ashley, 
Ritchie, and Brandt |Tlj . followed by quantum pertur- 
bation theory [T^ . and later further perturbative Born 
terms jlj, 13 have also been considered. Exact numer- 
ical solutions to the time-dependent quantum problem 
have been demonstrated by Mikkelsen and Flyvbjerg [l5j . 
The use of a harmonic oscillator model, as was previ- 
ously argued is justified for reasons of simplicity 
and at the same time good quality of results in compar- 
ison with observed stopping powers. The preset work 
adds yet another reason, as it demonstrates how classical 
results such as 0, 0] can be used to obtain exact quan- 
tum answers. Thus, the approach presented here does 
not involve heavy numerical calculations. It presents a 
simple and transparent way to understand and to treat 
exactly the Coulomb excitation of an oscillator. As an 
application of the presented technique, rather than re- 
peating previous numerical calculations, we concentrate 
on a less explored aspect of Coulomb excitation; here the 



overall excitation probability, the dipole and quadrupole 
transition amplitudes, and their interplay are discussed. 

We use a connection between classical and quantum 
mechanics which is usually complex and in most cases 
takes an approximate quasiclassical form. However, in 
situations where the dynamics of a system is deter- 
mined by linear equations of motion, with arbitrary time- 
dependent parameters, the classical-to-quantum corre- 
spondence is exact. The harmonic oscillator which is 
driven by the time-dependence of its parameters, such as 
mass, frequency or external force, is an example of such 
a situation. Most theoretical developments in the field 
of parametric excitations of harmonic systems were done 
more than three decades ago [H[ll El H . Despite 
this, it is only during the past decade theorists began to 
embrace the beauty of this problem and to discuss its rel- 
evance to many physical processes. Photon generation, 
lasers 22], quantum interference and two-photon excita- 
tions |23|. phase transitions, critical phenomena 0, HH j 
information theory |26| . ion traps, oriented chiral conden- 
sates [27j are examples of recent developments. Squeezed 
states that appear as solutions in parametrically driven 
harmonic systems, also known as two-photon coherent 
states, generalize the concept of coherent states I2II t o in- 
clude two-quanta couplings in the radiation field |29ll3fj| . 

By presenting this work we further expand applications 
of the squeezed state formalism to problems of Coulomb 
excitation. We focus our efforts on classical and quan- 
tum effects resulting from the interplay of one- (dipole) 
and two-quanta (quadrupole) couplings of an oscillator 
system to an external field. Using symmetries we derive 
simple formulas for excitation probabilities and transition 
amplitudes that explicitly show quantum interference ef- 
fects and can be used to improve perturbative calcula- 
tions. 

This work is divided into four sections. In Sec. ITT1 we 
review the properties of harmonic-oscillator systems and 
the transition from a classical to a quantum description, 
which is based on canonical transformation. We summa- 
rize the properties of the parameters that describe the 
initial state to final state transition (S'-matrix) and de- 
rive transition amplitudes and excitation probability rel- 
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evant to Coulomb excitation. Section Hill is dedicated to 
applications. Here we consider the classical solution and 
address the transition to the quantum treatment. We 
discuss quantum corrections, perturbative limits, charge 
asymmetry, and two-phonon excitations. The role of El 
and E2 interplay is emphasized in these discussions. We 
summarize our findings in Sec. IIVI 



II. PARAMETRIC EXCITATION OF COUPLED 
OSCILLATOR SYSTEM 

We consider a system of N oscillators with the unper- 
turbed Hamiltonian 
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(1) 



These oscillators are disturbed by a time-dependent force 
Fi (t) and time-dependent frequency parameters Qi i* (t) = 
Qv lit) , expressed by 



(2) 



It is assumed that flint (i) = in the infinite past t — > — oo 
and in the infinite future t — > oo . First we quantize this 
problem in the infinite past. Here 
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= ( ai e- iuj,t + aje^A 
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(3a) 
(3b) 



where otj and a; are time-independent creation and anni- 
hilation operators for the l-th mode. For any other time 
we introduce creation and annihilation operators b\(t) 
and bi(t) so that the coordinates and momenta in the 
Heisenberg representation are 



Pi 



( w "' 



— iuJit 



b\(t)e 
-bUt)e 



iujlt 
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The operators b'(t) and b(t) are defined in the interac- 
tion representation, as the time dependence due to flo is 
introduced explicitly in Eq. Q). For the infinite past we 
have by definition 

6(-oo) = a, 6 t (-oo)=a t . (5) 

At t — > +oo 6(i) and frt (£) also become time-independent 

6(oo) = 6, 6 t (oo)=6 t . (6) 

The transformation from the original operators a and 
to the intermediate b and w is determined by a general 
Bogoliubov transformation 



M*) = (**»' (*)«!' +«II'(*)4) + <*(*) 



(7) 



In order to simplify the notation we will write this in 
matrix form as 



b(t) = u(t)a + v(t)d* + c(t) 



(8) 



The perturbation H m t drives the nontrivial time de- 
pendence of these creation and annihilation operators 



i—b = [b(t), fl^] = T{t) + QT(t)b + e+(t)6+ , (9) 



where matrices 



U!l> 



(10) 



and a vector 



are introduced. The Q~ and Q + are hermitian and sym- 
metric matrices, respectively. The time dependence of 
the parameters u and v , as follows from @ , are deter- 
mined from 



and 



d 

i— u - 
dt 

d 

i— v 
dt 
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i— c - 
dt 



--Q-u + Q+v*, 
= Q-v + Q+ u* , 

Q-c+Q+c*+T 



(11a) 
(lib) 

(12) 



These equations are solved numerically with initial con- 
ditions 



u(— oo) = 1 , v(— oo) = , c(— oo) = , 



(13) 



which follow from JSJ. We stress here that solving Eqs. 
(|llfl and (|12|) is equivalent to finding a full classical solu- 
tion to the problem. This follows from the fact that the 
equations of motion are linear, so the time-dependence 
of a Heisenberg operator becomes identical to the time- 
dependence of its classical analog. 



A. Properties of the S-matrix 

In this work we will be interested only in the initial to 
final state transition probabilities, i.e., in the S matrix 
that comes from the transformation of operators 

b — ua + va' + c , b^ — u*^ + v*a + c* . (14) 

The inverse transformation is 

a = vtb - w T 6 f + c' , a f = u T tf - + c'* , (15) 

where 

(16) 



/ t i T * 

c = —u'c + v c 
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Transformation i|14|) is canonical and preserves commuta- 
tion relations so that if [a/, a\,] = Sw then {b[, b\,] = 5u> ; 
this puts some conditions on the matrices u and v 



uu' 
,t _ 



vv' 



T * 

uu' — V V 



uv — vu = U , 
t T * n 

u'v — v u = U . 



(17a) 
(17b) 



These linear canonical transformations, so-called sym- 
plectic transformations, are often encountered in classi- 
cal mechanics. The generators of the above symplectic 
transformation form a grou p which generalizes the sim- 
ple SU(1, 1) group |3ll|32,|33| 01 a one-dimensional os- 
cillator. In general any infinitesimal time evolution can 
be described by a linear canonical transformation which 
preserves the phase space in accordance with Liouville's 
theorem. The beauty of our case is that the linearity re- 
mains for any period of time evolution. This allows for 
an exact classical-to-quantum correspondence. Besides 
the method that we implement in this paper, the Wigner 
transformation [34| that establishes the correspondence 



between the classical phase space distribution and the 
quantum density matrix can also be used to exactly re- 
construct the quantum solution from its classical analog 
|35|. The application of the Wigner transformation in 
studies of various properties of the harmonic oscillator 
has been demonstrated repeatedly [2fil I3H IstI ] . 



B. Survival probability 

Here we assume that u , v and c are known. Using 
these we determine transition amplitudes between quan- 
tum states. Coherent states, 



l«)=e- W2/2 E^f 



(18) 



provide the best means for addressin g th is problem. In 
this basis the evolution operator U is [38| 



(/?|£/|a)=^exp(-i(|/? 2 | + |a 2 |) 



ia T (u*)- 1 wa + /3 t (u t )- 1 ( 



•/3 t (t* t )-V-a T (u*) 



*\-i c * 



r 



(19) 



The normalization constant A is determined by the com- 
pleteness condition 



d™PMU\a)\* = l, 



(20) 



where integration goes through both real and imaginary 
parts of (3 . This is an important parameter, since it sets 
the scale of all transitions and it determines the proba- 
bility to remain in the ground state Po,o = |(0|f/|0)| 2 = 
\A°\ . The quantity P , 

is given by the Gaussian inte- 
gral, where a is set to zero for simplicity 



7T iV f 

— = / d 2N (3exp {-|/3| 2 + 3£(/3 T p/3 + 2c T /3 - 2c t p/3)} 
"0,0 J 

(21) 



We introduce here the matrix 
p = v{u*)~ 



(22) 



which is a characteristic of quadrupole excitation. Be- 
sides some trivial normalization that will be further dis- 
cussed, the matrix element pi j is the amplitude of a single 
step two-quanta excitation of z-th and j-th modes. The 
classical time evolution of this matrix follows from the 
time evolution of u and v ljll|l . and is given by 



Q~P + pQ + p+Q^ 



(23) 



This result explicitly shows that p is a symmetric matrix 
(P = P T )- 



The multidimensional Gaussian integral can be taken 



d 2 "bexp(-b T Ab + 2s T b) 



where relevant to our case 



VdetA 



exp (s A 1 s) 



(24) 



, &/3 \ _ ( a?c A ( 1 - dip 3p 

~ 1 9/3 J ' C ~ I 3?c I ' A _ I 3p 1 + dip 



Q T c, 



Q = 



-dip 
-Qp 



xsp 

-(1 + ftp) 



(25) 



(26) 



Because of the previously discussed properties of the Bo- 
goliubov transformation, inversion of the matrix A is 
not needed. It can be shown that — c T QA : Q T c = 
-etc + di(c T p*c) . This results in 



Poo = V dct(l - ptp) exp [-c f c - 



(c T p*c)] . (27) 



The survival probability in the above form unifies the 
roles of El and E2 processes, and establishes a direct and 
simple relation between the quantities of perturbation 
theory, p and c , and the exact answer. There is a clear 
meaning of the different terms in the product (|27l) . The 
cxp(— c^c) represents the excitation of a coherent state 
due to the dipole interaction; the term y/det(l — p*p) 
is a result of quadrupole squeezing p7T |. The remaining 
part exp[5R(c T p* c)] is a result of quantum interference. 
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C. Transition Amplitudes 

Equation l|19f) can be used to obtain transition am- 
plitudes. Coefficients in a Taylor expansion of (|19|) are 
related to transition amplitudes between oscillator states 
|27|. with one-dimensional oscillator, for example, 



ifi\U\a) 



-CM 



-|/3| 2 )/2 Y,(n\U\m 



(28) 



The energy transfer, which is important for calculating 
the stopping power, can be calculated using Eq. (fHI) . For 
the excitation from the ground state we obtain 



A£=<0|5>^|0>=£ 



E 



Ul\VW 



(33) 



The above presents an average energy transfer. The en- 
ergy spread can also be calculated as in [27}. 



Below we denote by \i,j,k...) the normalized state con- 
taining single quanta in i,j,k... oscillator modes (later 
we will use directions x y and z). A single-quantum ex- 
citation amplitude is given by 



(j\U\0) = -A {( u t r y }j = yl°{(c — pc*)} 



3 > 



(29) 



where the subscript j denotes a component of the vector 
in the brackets. Since the interaction vanishes in both 
the infinite past and infinite future, the transition prob- 
abilities are symmetric with respect to initial and final 
states, in agreement with general principles of quantum 
mechanics. However, the decay amplitude from an ex- 
cited one-phonon state to the ground state, which is 



<0|J7|j) = -A°{(u*)-V} 



3 ' 



(30) 



differs from the excitation amplitude by a non trivial 
phase. This emphasizes the effect of quantum interfer- 
ence between dipole and quadrupole amplitudes. Setting 
either c or p to zero restores the symmetry. 

The transition amplitude from the ground state to the 
excited state with two oscillator quanta, requires expan- 
sion of (|19|l to the order 1 



D. One dimensional harmonic oscillator 

The special case of a one-dimensional harmonic oscilla- 
tor with time-dependent parameters has been extensively 
studied previously [H El E E El H E3 . According 
to Eq. l|2*7Jl the probability to remain in the ground state 
is 



^oo = Vl-IPl 2 cxp{-| 



where <A p and 



[l-|p|cos(^-20 c )]} , 
(34) 

„ are phases of p and c , respectively. This 
agrees with the previously obtained result ^3- We will 
further on encounter a one-dimensional case where c = . 
Here only even-quanta transitions are possible, and the 
excitation probabilities from the ground state and first 
excited state by 2n quanta are 



Po;. 



(2n) 



2 2 "(n!) : 



\p\*"y/l=\tf 



(35) 



Pi 



2n+l 



(2n + 1)! 
2 2 ™(n!) 2 



\p\ 2n (l-\p\ 2 ) 3/2 . (36) 



(ij\U\0) 



A° Pij + 



(i\u\o)W\Q) 

A° 



(31) 



The first term in the above equation is due to a direct 
quadrupole transition, since its amplitude is directly pro- 
portional to p . The second term describes the second- 
order dipole process. 

For an initially excited system the transition ampli- 
tudes can be calculated in a similar manner, by expand- 
ing in both a and (3 ; for example 



(z|c/|j) = aVX- 



A" 



(32) 



Tailor expansion of a Gaussian form can be used to 
define Hermite polynomials, in particular for the multi- 
variable case such as Eq. (IS}. This, in turn, allows for 
an analytic expression of transition amplitudes between 
oscillator states. Useful iterative relations between am- 
plitudes that follow from properties of Hermite polynomi- 
als can be obtained in this way, or alternatively, iterative 
relations can be derived directly from the properties of 
the evolution operator [3^| . 



III. COULOMB EXCITATION 

We consider here the problem of Coulomb excitation 
of a particle bound in a harmonic-oscillator potential, by 
a projectile of charge Z. We assume that the projec- 
tile moves along a fixed trajectory R(t) in the xy plane. 
In order to utilize the treatment discussed above we ex- 
pand the Coulomb potential up to quadratic terms in the 
particle coordinate q , namely, treating it in the dipole- 
quadrupole approximation 



C 



c 



V covl (t) = - 3 q-R , 2R5 



(3(<f • R) 



2 - q 2 R 2 



(37) 



where C = Ze 2 . The non-zero components of H mt de- 
fined in Eq. © are 



F x (t) 



CR X 
R 3 ' 



Fy(t) 



Q, 



C_ 
R? 



i 3-fiL — R 



CRy 



Jj5 ( 3R l 



(38) 



R 2 
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This leads to 



Qxy — 



lyx 



3C R x Ry 

p? 



Q = _£. 



(39) 



These expressions show that the mode in the direction 
perpendicular to the scattering plane is decoupled, and it 
is excited by the quadrupole term only. Further, we as- 
sume a straight line trajectory so that R(t) = (b, vt, 0) . 
This does not simplify the problem, since our formal- 
ism can be used for any trajectory. This assumption, 
however, allows us to compare with previous results 
[IH \12L [l3j • We introduce the following notations 



3i = 



C 



c 



2b 2 v y/LOiUjj 



and 



6 



v 



(40) 



(41) 



where i and j are indices referring to x , y and z direc- 
tions. Often we will assume lu x = lu v = lu z = lo ; in that 
case the corresponding indices are omitted. The param- 
eters and q can be viewed as dipole and quadrupole 
strengths, so that their ratio 1/(&v2cj) determines the 
applicability of the multipole expansion. 



Perturbative treatment 



<k»=-2iZ x K 1 {£ x ), c« = 2£ r ff (£ J ,) , cW=0. 

(44) 

We note that Cy is real, and c£ is imaginary. 

Using the squeezing parameter p can be calcu- 
lated in a similar Born-type perturbative manner. In the 
lowest order we obtain 



Pxy 



Pyy = %i°iyy£,l K o{ 2 iy) » 



4q^ y^xly (£x +Zy)Kx (£* + £y) , 

Pzz = 4iHzz€zKi(2£ z ) . 



(45a) 
(45b) 

(45c) 
(45d) 



The above results can also be obtained using time de- 
pendent perturbation theory. The first order excitation 
probability in the i-th direction is P} — \ci\ 2 . The prob- 
ability to excite quanta i,j in a single step quadrupole 
process is to the lowest order Pj^ = (1 + Sij)\pij\ 2 /A . 
These results can be compared with the exact transition 
amplitudes in Eqs. (J2U) and l(3T|) . 

In the sudden limit £ — ► the dipole contribution 
comes from the transverse direction: c x sa — 2id , while for 
quadrupole excitations \p xx \ ~ \pzz\ ~ 2q are dominant 
in this limit. The adiabatic limit £ — > oo is a textbook 
example of preservation of adiabatic invariances with ex- 
ponential precision |45j . 



Here we will discuss several approximations. The first 
one is the usual approximation relevant to multipole ex- 
pansion of the Coulomb field. It utilizes the smallness 
of the photon wave length as compared to the size of 
the system. The next two limits are the "sudden limit" 
where (<1 and the "adiabatic limit" £ 3> 1 . Finally, we 
will use the term Born approximation which corresponds 
to some order in the iterative solution of the classical 
equations Ijlll) and l|T2(l . Since the classical equations of 
motion allow here for an exact reconstruction of the wave 
function, the use of this terminology does not contradict 
its quantum notion. 

All of the results in this study, such as in Eqs. (|27fl. 
(|29|l . and are fully determined by the parameters 

c and p , which may be referred to as shift (dipole) and 
squeezing (quadrupole) parameters, respectively. Both of 
these parameters can be expanded using the above ap- 
proximations. Each order of Born approximation results 
in the next order correction in terms of the quadrupole 
strength 1)40 [I. Thus an expansion of c and p is a series 

C = c (n) (0q n_1 . P= E r(,i) (0q", ( 42 ) 

n=l,2,... n=l,2,... 

where n = 1,2 . . . is the corresponding Born order. The 
coefficients c^™^ and are £ dependent only. 

The lowest order long-range dipole approximation is a 
textbooks example |44|- Using Eq. l|12fl we obtain 

c«0c« - -t^Le^dt. (43) 
,/-oo V2w/ 



B. Z A corrections in classical solutions 

The long range dipole contribution is the main part of 
the cross section, thus any corrections to c are the most 
important. We denote as co(i) the inhomogeneous part 
of the solution to Eq. l(T2l 

c (t) =-i j T{t')dt' . (46) 

J — oo 

The Co = Co(oo) = dc^ 1 is the previously discussed first 
Born dipole amplitude. The second Born correction, 
which here simply means the second order approxima- 
tion to the classical solution of Eq. i|12|) is 

/oo 
[Q-c (t) + Q + c*(t)]dt. (47) 
-oo 

This integral can be studied numerically or perturba- 
tively with respect to £ . In Fig. Q]the dependence of the 
real and imaginary parts of c are shown as functions of 
£ . In the adiabatic limit, which corresponds to large £ all 
of these functions decay exponentially: ~ exp(— 2£) 
as dictated by the adiabatic nature |4^. Approaching 
the sudden limit all corrections go to zero as powers of 
£ . In the £ — > limit Qxi 1 '' is the dominant term, and the 
lowest order correction to it is Sci 2 '' w 37r£ 2 . This leads 
to the following approximation for the probability |ll| 

P x = 4c) 2 (l-37rq£ 2 ). (48) 
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FIG. 1: The lowest order polarization corrections to the dipole 
amplitude are plotted as a function of £ . The insert shows the 
behavior of Q(c^ ) and $l(c^ ) relative to 3-7r£ 2 . 

The above correction behaves as Z 3 with the charge of 
the projectile, and reduces the excitation probability in 
repulsive Coulomb scattering. This correction is com- 
pletely classical as it appears in a solution of classical 
differential equations, see Note that the homoge- 

neous part of HI 2(1 which is the source of this polarization 
effect, is caused by the E2 field. 

Higher order corrections to the squeezing parameter 
p are relatively large but are generally less important 
due to the overall smallness of the quadrupole amplitude. 
Here again we iteratively solve the classical equation ((23(1 . 
This equation contains no dipole parameters, thus the 
correction on p is self-induced and would not change in 
the absence of the El field. Assuming that 

/t poo 
Q+(t')dt'= / dt'Q(t-t')Q+(t'), (49) 
-OO J — OO 

the second Born correction is 

iSp = [ dt(po(t)Q-{t) + Q-*(t)po{tj) . (50) 

For the case of equal frequencies, this expression can be 
calculated analytically 

8 Pxx = Sq 2 ? [^(20 - 20K 2 (20] - ^q 2 ^(l - 20e~ 2 « ,(51a) 
5p xy = -27rq 2 e 2 e- 2 « - 8iq 2 £ 2 ^i (20 , (51b) 

6 Pyy = -i^q 2 £(l + 20e~ 2 « , (51c) 
16 

Sp„ = -8q 2 ^i (20 + — iq 2 eK! (20 . (51d) 

In the £ — » limit, excitations due to the quadrupole 
field comes from the x and z directions, leading to the 



following approximate E2 contribution to the excitation 
probability 

AP p = 2q 2 (l + 7rq£(l - 40) + 2q 2 + . (52) 

The charge asymmetry in this case is of opposite sign, 
resulting in an enhanced probability for repulsive inter- 
action. This enhancement can only be observed in mea- 
surements if the dominant dipole transitions are blocked 
or restricted. 



C. Exact treatment 

We start this section by showing the comparison of 
the exact excitation probability and the first Born dipole 
approximation, Pborn = c o c o j see Fig. [21 The harmonic- 
oscillator model is rather crude and can not be used to 
fully understand features of the realistic Coulomb dissoci- 
ation of nuclei. However, it can still provide a qualitative 
picture. For example the Coulomb dissociation of | 7 F 
can be looked at assuming that the proton outside the 
oxygen core is bound by a harmonic potential. The fre- 
quency of the oscillator, related to the dipole excitation 
energy, is assumed to be of the order of 1 MeV. In this 
case a unit of the oscillator length corresponds to 6.4 fm. 
Assuming Coulomb excitation by ^Ni we have C = 6.3 . 
The choice of parameters used in some figures below is 
guided by this model. In Fig. [21 we use impact parameter 
b = 4 (b — 25.6 fm), C = 6 , and the range of velocities 
from 0.2 to 20 (oscillator units) that corresponds to the 
incident energies from 0.02 to 200 MeV/u. The choice of 
impact parameter is restricted by the range of applicabil- 
ity of the multipole expansion \J~2b 3> 1 . Our discussions 
also imply that the probability to find the nucleon in the 
oscillator at a distance from origin larger than impact 
parameter is negligible. For 6 = 2 this probability is 4%, 
assuming the ground state wave function. However for 
b = 4 that we use, this is less than 10 -3 %. 

As seen in Fig. [21 the effects of higher order corrections 
are significant. The a difference between repulsive and 
attractive interactions is also transparent. This confirms 
the recent results of more realistic calculations [l| . 

1. Higher-order corrections 

From the classical solutions for p and c exact quantum 
results can be extracted using Eq. (|2"7|l . In order to 
discuss the impact of a quantum treatment we can make 
an expansion of the exact answer in Eq. ((27(1 , which gives 

P x = I- Poo » c^c-^- + ^TT(p^p)-^(c T p*c) (53) 
or ii + AP C + AP P + AP inter . 
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FIG. 2: Excitation probability as a function of incident ve- 
locity, at the impact parameter b — 4 . The dashed line is the 
result of the repulsive interaction, (7 = 6. The dotted line cor- 
responds to the attractive interaction C = —6 . These results 
are compared to the first Born dipole approximation shown 
by the solid line. 



FIG. 3: Corrections to the first Born El probability. Values 
are relative to Pbom- Logarithmic scale is used to show contri- 
butions of very different magnitude. Corrections AP c /Pborn 
AP P / Pbom and APmtcr/Pbom are shown in dotted, dashed 
and solid lines, respectively. Parameters 6 = 4 and C = 6 are 
used for this figure. 



We identify here three main corrections to the first Born 
dipole excitation probability. The first is 



AP r 



c'c - 



(etc) 2 



p, 



born ■ 



which is a correction to the dipole transition probability. 
It is known that the excitation of a quantum oscillator 
by an external force creates a coherent state. The dipole 
amplitude c , which is related to a classical shift of the 
oscillator from its equilibrium position, appears in the 
exponent of Eq. H27[) . This accounts for all second and 
higher order El processes. Exponentiation of the dipole 
amplitude, in order to account for the loss of strength to 
multi-phonon excitations, has been previously discussed 
[ir| . We stress here that other parts of AP C are of clas- 
sical origin and mainly Z 3 -dependent. They come from 
the effect of the quadrupole field on the shift of the os- 
cillator from its equilibrium. This creates a deviation of 
c from the Born result cq as discussed in Sec. IIII Bl 
The next term in (|53[) . 

AP P = Tr(pt p )/2 , 

is a quadrupole contribution, which comes from the ex- 
pansion of the square root in i|27|) . Similar to the ex- 
ponent in the dipole term, the use of the square root 
accounts for all higher order quadrupole transitions. 

Finally, a very interesting and purely quantum inter- 
ference effect appears through the last term in 15311 . 



Although it may seem that this is a quantum Barkas cor- 
rection, correction vanishes to lowest order. This follows 
from Eqs. lO and gSJ). 

In Fig. |3|the absolute values of AP C , AP p , and APi ntor 
relative to the first Born El probability are plotted as 
functions of velocity. The figure demonstrates that domi- 
nant corrections to the first order dipole probability come 
from second and higher order dipole transitions, and from 
the lowest order quadrupole. The interference term con- 
tributes only on the level of 1%, and quickly vanishes in 
the sudden limit. As discussed above, the interference 
term vanishes to order Z 3 . Thus it is a Z 4 term to lead- 
ing order. 



2. Z -dependence 

In order to examine the Z-odd part of the excitation 
probability from a different point of view we introduce 
the Barkas factor 



B 



PAZ) 



-z) 



P X (Z) + P x (-Z) 



(54) 



AR, 



c p c). 



As discussed above this effect owes its existence to the 
classical dynamic interplay between El and E2 fields. 
With the validity of a multipole expansion, the bulk 
of this correction is expected to be proportional to the 
quadrupole amplitude q. In Fig. ||B/q is plotted as a 
function of £ . Various cases corresponding to different 
impact parameters and projectile charges are considered. 
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FIG. 4: The Barkas factor as a function of £ , for various pa- 
rameters. The classical correction to the sudden limit, 37r£ 2 , 
is shown by a dashed curve. 



FIG. 5: The transition probability from the ground state to 
the two-quantum excitated state in the z direction is shown 
as a function of incident velocity. Solid and dotted lines corre- 
spond to repulsive and attractive Coulomb force, respectively. 



In the second Born approximation, the B/q is only a 
function of £ , thus independent of both charge and im- 
pact parameter [ljj. Deviations observed on the plot 
emphasize the significance of other corrections that scale 
as higher powers in q and higher odd powers in Z . 

The charge asymmetry related to the dipole- 
quadrupole interference appears naturally in various 
transition amplitudes. The Fig. [S] shows the transi- 
tion probability from the ground state to the state with 
two quanta in z direction. In this transition the vanish- 
ing out-of-plane dipole force makes the quadrupole term 
dominant. As demonstrated by Eq. i|52|) the Z 3 Barkas- 
type charge dependence of the quadrupole term p is of 
the opposite sign in comparison to the dipole-quadrupole 
case shown in Figs. |2land0] In agreement with this, an 
enhancement of the excitation probability in repulsive 
Coulomb interaction is observed in Fig. [3] 



3. Interplay of dipole and quadrupole fields 

Our harmonic-oscillator model incorporates El and E2 
electromagnetic processes exactly. The role and interplay 
of these processes in multi-phonon excitation may be of 
interest in the physics of Giant Dipole Resonances (GDR) 
[4^. |47| . Although precise calculations with inclusion of 
the relevant microscopic physics of GDR are beyond the 
scope of this paper, we consider here the effect that the 
E2 process plays in the two phonon excitation probability 



We introduce a quadrupole coupling strength x which 
simply scales the quadrupole part of the perturbation 
Hamiltonian J2J. This can be viewed as scaling the 
quadrupole strength q — ► xq . The limits x — and 
x = 1 naturally correspond to the cases when E2 is ne- 
glected and when E2 is at its natural value, respectively. 
Eq. as one would expect from perturbation the- 

ory, shows that the amplitude for a two-phonon excita- 
tion is given by the sum of two terms that reflect the 
first order quadrupole and second order dipole processes. 
The El and E2 interplay is not trivial here and comes 
in several places. The first term in (|31fl at large quadru- 
ple strength is directly proportional to x . Thus in this 
limit a quadratic scaling of the excitation amplitude is 
expected. At a smaller strength a substantial contribu- 
tion to the second term in i|31[l appears as a classical 
dynamic polarization effect, influencing the amplitude c . 
This can reduce the excitation probability for a repulsive 
interaction. Finally, the phases of terms in 131|) can in- 
terfere. It must be noted that the dynamic polarization 
disappears at very low and very high velocities. Also 
both X) ~ 1/v and q ~ 1/v but since dipole comes in 
the second order, the quadrupole amplitude dominates 
at high velocity leading to 

P ,2(z)=4qV. (55) 

In Fig. EJthe two-phonon excitation probability as a func- 
tion of the coupling strength x is shown relative to the 
case with no E2, x = . Observed features are in agree- 
ment with the above discussion. Cases corresponding to 
repulsive interaction, shown by solid lines, always lead to 
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FIG. 6: Relative two-phonon excitation probability as a func- 
tion of the quadrupole strength. 

reduced probability as compared to kinematically identi- 
cal situations but attractive interaction (dashed lines). 
The dynamic polarization effect is maximized around 
v = 3 which is consistent with the enhanced charge asym- 
metry observed in Fig. In this velocity regime, the 
polarization effect is so strong that for repulsive interac- 
tions "turning on" the quadrupole field (x — I) leads to 
a reduction of a two-phonon excitation probability. 

IV. CONCLUSIONS 

In this work we considered a simple model of a charged 
particle in a harmonic-oscillator potential. The system 
is excited by the electric dipole and quadrupole fields 
of a charged projectile passing by. We have devel- 
oped an application of the squeezed-state formalism to 
solve this problem exactly. In terms of the one-phonon 
dipole excitation amplitudes c and the direct two-phonon 
quadrupole excitation amplitudes p we derived simple 
expressions for various transition amplitudes, and the 
total and multi-phonon excitation probabilities. The pa- 



rameters c and p are determined from classical equations. 
Our results are exact for a harmonic oscillator. For a 
general case they can be used as a parameterization, ex- 
tending the higher-order (multi-phonon) dipole approxi- 
mation introduced by Norbury and Baur [46| to include 
the quadrupole field and related interference effects. 

We discussed the interplay between the electric dipole 
and quadrupole fields and the resulting effects on vari- 
ous excitation probabilities. In the intermediate range 
of energies, between instant and adiabatic limits, (0.5 < 
£ < 7) the influence of dynamic polarization dominates. 
It was shown that the origin of this effect is purely clas- 
sical. The quadrupole polarization influences the effect 
of the dipole force, leading to a significant change in the 
excitation probability. This phenomenon, known in stop- 
ping power theory as Barkas effect, in the lowest order 
leads to a Z z charge dependence, and reduces the excita- 
tion probability for repulsive Coulomb interactions. We 
found a similar classical effect of self-induced polarization 
in the quadrupole amplitude. This effect is of the oppo- 
site sign and for the same repulsive kinematics results in 
an enhancement of the quadrupole amplitude. The effect 
was demonstrated in z-polarized two-phonon excitations. 
Our exact results exhibit some additional higher-order 
quantum contributions, but corrections from these are 
small. We presented a number of numerical calculations 
and plots that demonstrate a variety of observable phe- 
nomena that can be attributed to the E1-E2 interplay. 

By presenting this work we hope to introduce the co- 
herent plus squeezed-state formalism to the field of nu- 
clear reaction physics and nuclear excitations. Numerous 
developments that use dipole excitations can be extended 
to include quadrupole transitions. Not every physical sit- 
uation can be described by a harmonic oscillator. How- 
ever, armed with an exact solution and a full set of time- 
dependent wave functions, perturbation theory based on 
squeezed states can be a promising future direction. 
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